Quantum theory of optical temporal phase and instantaneous frequency 
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We propose a general quantum theory ol optical phase and instantaneous frequency in the time 
domain for slowly varying optical signals. Guided by classical estimation theory, we design homodyne 
phase-locked loops that enable quantum-limited measurements of temporal phase and instantaneous 
frequency. Standard and Heisenberg quantum limits to such measurements are then derived. For 
optical sensing applications, we propose multipass and Fabry-Perot position and velocity sensors 
that take advantage of the signal-to-noise-ratio enhancement effect of wideband angle modulation 
without requiring nonclassical light. We also generalize our theory to three spatial dimensions for 
nonrelativistic bosons and define an Hermitian fluid velocity operator, which provides a theoretical 
underpinning to the current-algebra approach of quantum hydrodynamics. 
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I. INTRODUCTION 

Frequency modulation (FM) radio sounds better than 
amplitude modulation (AM) radio, because by encoding 
a message in the instantaneous frequency, defined as the 
rate of change of the signal phase, FM uses more trans- 
mission bandwidth to improve the signal-to-noise ratio 
(SNR) [l|, 0j 3- Ii^ optical communications, FM and 
phase modulation (PM) techniques have also received 
significant attention and can offer higher SNRs in pho- 
tonic links than standard intensity modulation [j . Quan- 
tum statistics is expected to play a major role in the 
performance of coherent optical communication systems 
@) S 0] ' but although a treatment of digital FM quantum 
noise has been briefiy mentioned by Yuen i7l| , a more fun- 
damental quantum theory of temporal phase and instan- 
taneous frequency is not yet available. In optical sensing, 
on the other hand, a considerable amount of research has 
been devoted to the study of ultimate quantum limits 
to phase [3, [3, (Tfl, [lH and velocity [lj| measurements. 
A quantum description of temporal phase and instanta- 
neous frequency is again needed to deal with any rapid 
change in the measured parameters. For example, in laser 
Doppler velocimetry, which has been widely used in fiuid 
dynamics [iJI and blood fiow diagnostics [ij, the veloc- 
ity of the interrogated sample actually shifts the instan- 
taneous frequency of the refiected optical signal and not 
the average frequency directly, so the qua ntum limits to 
velocity measurements derived in Ref. [I4I are no longer 
accurate when the velocity changes rapidly. 

A rigorous quantum description of the optical phase is, 
unfortunately, not trivial even for one optical mode, and 
has been a subject of long-running debate @, Ha M, Il3j 
[iSl . [ig . [20 . [2H since the issue was first raised by Dirac at 
the birth of quantum electrodynamics |22l | . The difficulty 
of defining a phase operator is shared by the superfluid 



and Bose-Einstein condensate community [23| , where the 
fluid velocity in Landau's formulation of quantum hydro- 
dynamics [2J] is often related to the phase gradient of an 
order parameter. Beyond the Bogoliubov approximation, 
while Landau has proposed a fluid velocity operator [2J] , 
many have argued against its existence [25| . 

In this paper, we show how the temporal phase and 
the instantaneous frequency of a slowly varying opti- 
cal signal can be treated consistently in the quantum 
regime. By limiting the bandwidth of the optical Hilbert 
space, we discretize the continuous time domain into or- 
thogonal "wave-packet" modes via the sampling theorem. 
This crucial step allows us to apply previous studies of 
discrete-mode quantum phase to the time domain. To 
enable quantum-limited temporal-phase measurements 
in practice, we use estimation theory to design homo- 
dyne phase-locked-loop measurement schemes, similar to 
adaptive phase measurements |27i, 2&, 2^ ■ Standard and 
Heisenberg quantum limits to the accuracy of temporal- 
phase and instantaneous-frequency measurements for co- 
herent states and squeezed states are then derived. For 
sensing applications, we propose multipass and Fabry- 
Perot position and velocity sensors that take advantage 
of the SNR enhancement effect of wideband angle modu- 
lation. Finally, we generalize our time-domain formal- 
ism for optics to three spatial dimensions for nonrel- 
ativistic bosons and define an Hermitian fluid velocity 
operator, which provides a theoretical underpinning for 
the current -algebra approach to quantum hydrodynamics 
[23,l2l[2i,l2i. 
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II. QUANTIZATION OF BANDLIMITED 
OPTICAL FIELDS 



To describe a typical free-space or fiber optical experi- 
ment, propagating modes are usually quantized in terms 
of the frequency, transverse momentum, and polarization 
degrees of freedom [^, |30, [Sll ■ We shall consider a plane 
wave of a certain polarization in free space or one fiber 



mode at an observation plane along the optical axis, and 
allow only the frequency degree of freedom for simplicity. 
In the slowly varying envelope regime, one can regard the 
optical fields as a one-dimensional nonrelativistic many- 
boson system and quantize accordingly [y, l30l . l3ll . |32| |. 
The positive-frequency electric-field operator can then be 
approximated as 



E'-+\t) ex i(t)exp(-i27r/ot), 



(2.1) 



where A{t) is the envelope annihilation operator and /p 
is the carrier optical frequency. We shall make the phys- 
ically reasonable assumption that our ability to create, 
manipulate, transmit, and measure the optical fields is 
limited to a certain bandwidth B near the carrier optical 
frequency /q. We can then express A(t) as the bandlim- 
ited Fourier transform of the frequency-domain annihila- 
tion operator a(/), 



B/2 

A{t)= I dfa{f)exp{~t2TTft). 



(2.2) 



We make the important observation that the time- 
domain commutator given by Eq. (|2.4p goes to zero at 
discrete intervals, as shown in Fig. (TJ when t — i' is a 
nonzero-integer multiple of l/B. If we discretize time 
into such intervals, that is, let 
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the commutator at these sampled times becomes 



[A{t,),AHh)]^BS,,^j^. 



(2.7) 



(2., 



The sampled times have become orthogonal "wave- 
packet" modes, as described qualitatively but not sub- 
stantiated in Refs. Q. The operators can be renormal- 
ized to give the conventional discrete-mode commutator, 



In other words, we only need the set of operators a(/) 
and a'l'(/) for which / falls within the band (|/| < B/2) 
to describe our experiments. To satisfy the slowly vary- 
ing envelope approximation, we require B/2 <gC /p. The 
commutation relation between a{f) and the correspond- 
ing creation operator a^f) is 



[a(/),at(/')]=<5(/-/'). 



(2.3) 



which leads to the following time-domain commutation 
relation, 

[A{t),A^{t')] ^BsmcB{t~t'), (2.4) 

where the sine function is defined as 



, , sin(7ra;) 
smc(a;) = 



TTX 



(2.5) 



aj = A[tj)\'St, 



[aj,al] =6 



'jk, 



(2.9) 



and the continuous-time and continuous-frequency oper- 
ators can be reconstructed using the sampling theorem. 



A{t) = Vb ^ aj sine B{t -tj), 



(2.10) 



1 °° B 

°(^) = ^W E ^^ exp(*27r/i,) for |/| < -. (2.11) 



By virtue of the sampling theorem, the use of the 
discrete-time operators on the bandlimited Hilbert space 
is completely equivalent to the use of continuous-time or 
continuous-frequency operators. 

The bandlimited Hilbert space can be spanned by 
photon-number states in discrete wave-packet modes, 
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FIG. 1: (Color online). sincB(i-i') 
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(2.12) 
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where Uj is the number of photons in the wave-packet 
mode at time tj, and for convenience we use boldface 
n to denote the vector {. . . , n^, rij+i, . . . }. Appendix |^ 
lists other notations that we shall use to describe algebra 
in the discrete time domain. For a pure state |^), the 
photon-number representation is 



C[n]^(n|vl/), E|C[n]p = l. 



(2.14) 



Another way of spanning the Hilbert space is via 
nonorthogonal multiphoton time-measurement eigen- 



states [la ISll la 

In, ■ ■ ■ :Tm 



)^^it(ri)...it(^^)|0), 



(2.15) 



i = ^ / dri . . .dTAr|Ti, . . . ,rAr)(Ti, . . .,Tn\. 
N=0'' 

(2.16) 
The multiphoton wave function is then 

V'iv(ri,...,Tjv)= (ti,...,tjv|*), (2.17) 

^ /dTi...dTAr|V'w(ri,...,Tjv)|' = l, (2.18) 



Af=0' 



which is related to the photon-number representation by 

M 






(2.19) 
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(2.20) 



where Piifj) is an element of a permutation of 
{ri, . . . ,rAr} and the sum over I is over all N\ possible 
permutations so that <^n is symmetric. The inverse rela- 
tion is 



C[n] 



N 



1/2 



N\5t 
X i}}n{- . . ,tj, . . . ,tj,tj+i, . . . ,ij+i, . . .), (2.21) 

nj terms nj + i terms 

where N = 'Y^-nj. For example, a coherent state can be 
defined in the continuous time domain as 

N 



\A{t)) = exp 



dtA{t)A^{t) 



(2.22) 



A{t)\A{t))^A{t)\A{t)), N: 



dt\A{t)\^, (2.23) 



where A{t) is a bandlimitcd envelope function and N 
is the average photon number. In terms of the wave- 
function representation, 



^AT! 



Using Eq. (j2.2ip , we obtain the photon- number represen- 
tation for a coherent state, 



a J =Aitj)VSt, (2.25) 



which shows that any coherent state can be written as a 
collection of independent coherent states in wave-packet 
modes. 



III. CANONICAL TEMPORAL-PHASE 

MEASUREMENTS 

For each wave-packet mode at time ij, we can define 
the nonunitary Susskind-Glogower exponential-phase op- 
erator [15] as 



£,^ 



a J a] 



A{t/)At{t,) 



-.A{t,). 



(3.1) 



Since £j depends only on the field operators at time tj, 
it can be regarded as an instantaneous phase operator. 
Despite the nonunitary nature of £j , its eigenstates form 
a nonorthogonal basis of the Hilbert space, 

|0)=^exp(m-</))|n>, 0o < ^ < 0o + 27r, (3.2) 



i= /i?0|0)(0|. 



z.0.n^ 



(3.3) 



The Susskind-Glogower eigenstates can thus be used to 
define a temporal-phase positive operator-valued mea- 
sure (POVM), also called a probability operator measure 
(POM) [i,|2S[2li, 



!![</,] = !</,) (01, j Dct>n[ct>\^l 



(3.4) 



which corresponds to a measurement of instantan- 
taneous optical phases {. . . , ^j,(/)-,_|_i, . . . } at times 
{. . . ,tj, tj+i, . . . }. It can be shown, by generalizing the 
single-mode treatment in Ref. 21], that the POVM given 
by Eq. (|3.4p corresponds to the optimal temporal-phase 
measurements for any periodic cost function with non- 
negative Fourier coefficients. Following Leonhardt et al. 
[18J . we shall refer to the phase POVM measurement 
as the canonical phase measurement. In the single- 
mode case, an experimental canonical phase measure- 
ment scheme has been proposed by Pregnell and Pegg 
[33| . While generalization of their scheme to the time 
domain is conceivable, it is beyond the scope of this pa- 
per to investigate experimental canonical temporal-phase 
measurement schemes in detail. 

The canonical temporal-phase probability density is 



p[</.] = Tr{pn[</,]}, j DM4>] = 1- 



(3.5) 



For a pure state, one can define a temporal-phase repre- 
sentation as the discrete Fourier transform of the photon- 
number representation. 



C[4>] = (01*) = ^exp(-zn • 0)C[n], 



(3.6) 



the magnitude squared of which gives the probabiHty 
density p(</>). 

One can also extend the Pegg-Barnett formahsm to the 
time domain by considering a Hilbert space spanned by 
finite- photon- number states \n) with < n < s. The 



commutator for Uj and a, becomes 



[aj,al] = S.jk[l- isj + l)\sj)j{. 



(3.7) 



although applying the commutator to physical states 



the 



limit recovers the usual commutator 



[aj,aj,]p = Sjk- A unitary Pegg-Barnett exponential- 
phase operator can be defined in the photon-number ba- 



exp(i0j) = Y^ \nj - l)j{nj\ + exp[i{sj + l)(j)oj]\sj)j(0\. 

rij — l 

(3.8) 

By taking the limit Sj ^ oo at the end of calculations, the 
Pegg-Barnett theory predicts the same phase statistics 
as the canonical phase statistics governed by Eq. (|3.5p . 
and the two theories can be regarded as equivalent and 
complementary [9|. 

To define an instantaneous-frequency operator, con- 
sider the classical definition 



Fit) ^ -^^m- 

2tt at 



(3.9) 



(j)(t) can be multivalued or even undefined when the in- 
tensity is zero, so to avoid this ambiguity we write the 
instantaneous frequency in terms of exp[i0(i)]. 



F{t) = ^|exp[z0(i)]^exphz0(i)] -c.c. 






- sin[0(i) - ^{t')] 



(3.10) 



where c.c. denotes complex conjugate. In the quantum 
regime, we can discretize F{t) and express it exactly in 
terms of the Pegg-Barnett operator, 



F =^ 

^ 2TT5t 



Ed. 



-k sm( 



M, (3.11) 



o) 



1 



ex:p{i4>j — i(j)k) — H.c. 



(3.12) 



H.c. denotes Hermitian conjugate, and dj-fc is called 
the differentiator 34] , the discrete impulse response that 
corresponds to differentiation. 



. _ r(-i)^-V(j--fc), j^fc, 



0, 



j = k. 



(3.13) 



Fj is Hermitian and well defined in the finite-photon- 
number bandlimited Hilbert space, but given the 
difficulty of canonical phase measurements in prac- 
tice, it might be even more difficult to measure the 
instantaneous-frequency operator, unless approximations 
are made. 



IV. MAXIMUM A POSTERIORI ESTIMATION 

It is difhcult to perform canonical phase measure- 
ments, and first-order field measurements, such as ho- 
modyne detection and heterodyne detection, are often 
preferred in practice. Statistics of heterodyne detection 
and its variants [l7| are governed by the Q distribution 
[isl . [35L [so, [331 , which is broader than the Wigner dis- 
tribution that governs homodyne detection [3^, |3y|, so 
ideally one would like to use only homodyne detection 
to measure the phase and instantaneous frequency. In 
this section we would like to show how this can be done 
from the perspective of estimation theory, which natu- 
rally leads to the use of homodyne phase-locked loops to 
perform angle demodulation. 

The Wigner distribution is defined as [3y| 

W[x,y]= I Dbexp{b* ■a-ba*)x[b], (4.1) 

ncPbi \ , , , , 

^, a=-{x^iy), 4.2 

TT^ 2 

i 
-[b] = Tr|pexp('b-at-b*.a')|, (4.3) 



XI 



which governs the statistics of operators that can be ex- 
pressed in Weyl-ordered quadrature operators [1^, [30| , 

(/W[&,y]) = J DxDyf[x,y]W[x,y], (4.4) 

X = a + a/ , y = —i{a — a*), (4.5) 

DxDy = Y\ dxjdy.j, (4.6) 



where the superscript iW) denotes Weyl ordering. In 
the following we shall consider only quantum states that 
have nonnegative Wigner distributions, and assume that 
measurements of Weyl-ordered quadrature operators can 
always be realized. The Wigner distribution is then a 
qualified classical probability distribution that we can use 
in classical estimation theory. 

If the range of phase modulation exceeds 27r, one must 
unwrap the measured phase to unambiguously decode 
the message contained within. Phase unwrapping can 
be done only if one has a priori information about the 
message, because otherwise one would have no other 
way to distinguish <j)j from (f)j + 2n7r, n being an arbi- 
trary integer, in a measurement. To incorporate a priori 
information in angle demodulation, we shall use maxi- 
mum a posteriori (MAP) estimation [l|, [2]- Given an 
a priori probability distribution of the message P[m] — 
P{. ■ ■ , TTij, rrij+i, . . . ), MAP estimation seeks the message 
that maximizes the a posteriori probability distribution 



P[m\x,y] 



W[x,y\m]P[m] 
P[x,y] 



(4.7) 



for a set of measurements in terms of x and y, and is 
asymptotically efficient. Taking the logarithm of the a 



posteriori distribution and differentiating with respect to 
m, we arrive at the vectorial MAP equation, the solution 
of which gives the MAP estimate m, 



\7mlnW[x,y\m]+\7mlnP[m]) _=0, (4.8) 

^ / m—m 

where the gradient operator Vm, is defined as 

d d 



drrij ' 9toj_|_i ' 



(4.9) 



The MAP equation (|4.8p can be significantly simplified if 
the probability distributions are Gaussian. We therefore 
model the message m as a zero-mean Gaussian random 
process, 



/I 
P[m\ ex exp m ■ K^ 

Km = (Am® Am), 



m 



(4.10) 
(4.11) 



where Krn is the message covariance matrix. In angle 
modulation systems, the mean phase is a linear transfor- 
mation of the message. 



4> ^ H ■ m. 



(4.12) 



where H is a real impulse response. For example, the 
impulse response for PM is 



H^pl, 



(4.13) 



where /3 is called the modulation index and I is the iden- 
tity matrix, and the impulse response for FM is 



Hjk = -27r.F 



dt sine B{t — ife). 



(4.14) 



where J^ is called the frequency deviation. Although the 
Wigner distribution is Gaussian for squeezed states [Sq . 
we shall first study the simpler coherent states. 



A. Coherent states 

To derive the Wigner distribution for a coherent state, 
let us start with that for the vacuum, 



1 1 

Wo[xo,ya] ex exp ( --a;o • xq - -y^ ■ y^ 



(4.15) 



A coherent state with mean phase (f) and constant mean 
amplitude |a| is obtained by displacing the vacuum along 
the a^o quadrature followed by phase modulation. The 
resulting Wigner distribution is 

W[x,y\4)\ = Wn[xo[x,y],y^^[x,y]], (4.16) 

a;o = a;cos0 -I- y sin0 — 2|a|, (4-17) 

j/o — —X sin (/) + y cos (j), (4-18) 



which can be used as the conditional distribution in the 
MAP equation (|4.8|) . After some algebra, the MAP equa- 
tion becomes 



m = 2\a\Km-H' p, 



(4.19) 



where p is a quadrature field that depends linearly on the 
complex field a and nonlinearly on the estimated message 
m, 

p = —i[a cxp{—i(f)) — a* exp{i<p)] , (p = H ■ rh. 

(4.20) 
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FIG. 2: (Color online). A block diagram that represents the 
MAP estimation process. D{\a\) denotes displacement of the 
vacuum along the xq quadrature, and expicj) denotes phase 
modulation. 



The nonlinear MAP equation (|4.19|) can be represented 
conceptually by a block diagram, as shown in Fig. [21 
but cannot be solved analytically in general. Although 
the block diagram seems to suggest that homodyne mea- 
surements of the quadrature p in a feedback system can 
be used to produce the estimated message m, p de- 
pends recursively on </), and the impulse response matrix 
H ■ Km ■ H that transforms p to (f) is not necessarily 
causal. Even if one attempts to solve the MAP equation 
numerically, one needs the values of both quadratures of 
a at all times, so one would have to perform heterodyne 
detection and start with the Q distribution instead of the 
Wigner distribution. 

To solve Eq. (|4.19p approximately, let us linearize it so 
that the right-hand side of Eq. (|4.19p depends linearly on 
m and we can use well established linear estimation tech- 
niques to perform MAP estimation. Writing the complex 
field a in terms of Xq and yg ^^ 



a — exp(i0) \a 



Xq + ly^ 



(4.21) 



p can be rewritten as 

p = 2|a|sin(0-0) + 2;o, (4.22) 

where we have defined another quadrature field zq as 

2;o = a;osin(^-0) +yoCos(^-0). (4.23) 

Quantum noise enters the MAP estimation through the 
2;o quadrature. Because the vacuum Wigner distribution 
given by Eq. (j4.15p is invariant with respect to rotations 
in phase space, any quadrature in terms of Xq and y^ has 
the same Gaussian statistics as Xq or y^ and its statistics 



are independent of cf) or </>. Zg can therefore be regarded 
as an independent noise term. To linearize p, we assume 
that the estimated phase is at all times close to the mean 
phase of the field in the mean-square sense, 

((^-0)2)«1, (4.24) 

so that we can make the first-order approximation 

p^2\a\{^-4>) + Zo. (4.25) 

A linearized MAP equation can thus be obtained, 



B. Phase-locked loops as angle demodulators 



m«4|arK^ -iJ • U>-^ 



2bl 



(4.26) 



In this linear regime the MAP estimation is efficient. Let 
us define 



- Zq ~ 

2 a 



(4.27) 



so that the linearized MAP equation (|4.26p can be writ- 
ten as a pair of equations. 



m = A\a\'^Km.- H'^ • e, 



(4.28) 
(j)-Hm = e. (4.29) 

Multiplying Eq. (|i:^ by H and adding it to Eq. (g^H]), 



0= {4\a\^H-Km-H^ + 1) -e. 



(4.30) 



Comparing Eq. (|430l) with Eq. (|4?28l) . we find that m 
can be written as a linear transformation of (p, 

Tn = G(j}, (4.31) 

where G is the solution of the following equation, 

G ■ (4|apff ■Krr.-H^ + I)^ Ma\^Km ■ H^ ■ (4.32) 

To solve this equation, we assume that the message co- 
variance matrix Km is stationary, so that we can write 

{Km)jk=Krr,^0-k. (4.33) 

and define its power spectral density as 

SmU) = Y. ^-.J exp(z27r/i,). (4.34) 

3 

If H is time-invariant, 

Hjk = ffj-fc, H{f) = Y, H, expii2TTftj), (4.35) 



we can solve Eq. (I4.32p by Fourier transform and obtain 
a time-invariant solution for G, 



G{f) = 



Gj-k 



4\a\^Srr.{f)H*{f) 

4|a|25^(/)|iJ(/)|2 + l' 

1 /-^/^ 

B 



(4.36) 



dfG{f)exp[-i27rf{t,-tk)]. (4.37) 



-B/2 



Gj-k is called the optimum filter in linear estimation 
theory [i|, [2]- In the next section we shall show how 
the optimum filter can be implemented by a homodyne 
phase- locked loop. 
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FIG. 3: (Color online). A homodyne phase- locked loop. LO 
denotes local oscillator. 

Consider the homodyne phase-locked loop shown in 
Fig. [31 The output of the homodyne detection is given 
by 

p' = —i[aexp{—icf}') — a* cxp(i0')] (4.38) 

= 2|a|sin(0-</)') + z', (4.39) 

where z' is another quadrature field defined as 

z' = Xq sin(0 - (j)') + yo cos(0 - 0'), (4.40) 

which, like Zq, has Gaussian statistics independent of cf) 
and (p . If we again linearize p' by making the assump- 
tion, to be justified later, that the local-oscillator phase 
(f> is at all times close to the mean phase of the incoming 
field 4>, 



((^-0')'>«1, 
p' is given by 

p'^2\a\{4>-cP'), = 0- 



z 

2W 



0' is related to p' by the causal loop filter L, 

4>' = L-p'^L-2\a\{ct>-(j)'), 



(4.41) 



(4.42) 



(4.43) 



and (p can be written in terms of cp and a closed-loop 
filter L', 

4,' TuL' ■ 0, (4.44) 

where the closed-loop filter satisfy the equation 

{I + 2\a\L)-L' = 2\a\L. (4.45) 

To minimize {{cp — (p')'^) and ensure that Eq. (|4.4ip holds, 
we demand (p to be the minimum-mean-square-error es- 
timate of 0, using only past and present values of p' to 
ensure causality. For stationary 0, 0, and z' , L' is sim- 
ply the Wiener filter and can be solved by a well-known 
frequency-domain technique [l|, ylj as briefly described 
below. L'.f. is the solution of the discrete Wiener-Hopf 



equation, 

oo 

YL'kU,-k^V,, 



k=0 



U = A\afH-K 
V 



H' 



A\a\^H-Km-H^. 



(4.46) 

(4.47) 
(4.48) 



The Fourier transform of Uj-k can be factored into the 
form 



U{f)=X{f)X*{f), 



(4.49) 



where X{f) and l/X{f) are causal filters, if U{f) is a 
rational spectral density. Defining 



Vjf) 

[x*{f) 



1 

B 



-B/2 ^ U ) j.^o 

(4.50) 



the Fourier transform of i'_j, is given by 



L'if) 



1 



X{f) 
and the loop filter is then 

L{f) 



Vjf) 
X*{f) 



L'U) 



(4.51) 



(4.52) 



2\a\[l-L'{f)Y 

It can be shown that the loop filter L is also causal given a 
causal L' [l|, 13 • After passing through the phase- locked 
loop and the post-loop filter L , the output estimated 
message becomes 



L' -cj). 



(4.53) 



This equation is the same as the linearized MAP equation 
glsTjl if 

G = L" L'. (4.54) 

The post-loop filter in the frequency domain is thus given 

by 



L"{f) 



Gjf) 
L'ifY 



(4.55) 



To make the post-loop filter realizable, we allow delay in 
the estimated message, that is, we let 



h = Y.H'^- 



kirik+d, (f>: 



= E^.- 



kUlk+d, 



(4.56) 



and L can be made causal in the limit of infinite delay 
d [llj Q- III practice dSt only needs to be a few times 
larger than the correlation time of the message for L" 
to be well approximated by causal filters. Hence, the 
homodyne phase-locked loop is able to realize the MAP 
estimation for coherent states, provided that the approxi- 
mation given by Eq. (|4.4ip is valid and the local-oscillator 
phase follows closely the mean phase of the incoming field 
at all times. This result is hardly surprising, because a 
coherent state can be regarded as a classical signal with 
additive white Gaussian noise, in which case it is well 
known that the homodyne phase-locked loop realizes the 
optimum angle demodulator [l|, |2| . Our analysis so far 
closely mirrors that in classical estimation theory, except 
that we have worked with discrete time for consistency 
with previous sections. It is straightforward to take the 
continuous limit of our analysis in this section, if one so 
desires, by redefining some of the quantities and taking 
the limit B —^ oo. 



C. Squeezed states 

The MAP estimation analysis shows that quantum 
noise enters the estimation through the Zq quadrature, 
which depends on both Xq and j/q. If we squeeze the Zq 
quadrature of the vacuum before displacement and phase 
modulation, the noise entering the measurement of p can 
be reduced. The squeezed-vacuum Wigner distribution 
is 



Wo[xo,yf)\ cxexp 



2^0 ^1 -Co - 2^° ^^ ' ^° 



(4.57) 

(4.58) 

(4.59) 
(4.60) 



Co = xo cos{4> -(t))-yf^ sin(0 - 0), 
Zq = Xq sin(0 - 4>) + y„ cos(0 - 4>), 

Ki=(Co®Co>, K2={zo<»Zo). 

The MAP equation becomes 

m = 2\a\Km ■ H^ ■ K^^ ■ p, (4.61) 

as shown schematically in Fig. [H The linearized MAP 
equation is the same as Eq. (|4.3ip . but G now also de- 
pends on K2, 



G-{4\aYH-Kr 



H^ + K 



(4.62) 
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FIG. 4: (Color online). MAP estimation for squeezed states. 
S{zo) denotes squeezing of the zq quadrature. 

A homodyne phase-locked loop, with redesigned loop 
and post-loop filters L and L according to the G given 
by Eq. (|4.62p . can again be used to realize this MAP 
estimation, as shown in Fig. [51 although one must feed 
the local-oscillator phase <p' back to the quantum-state 
preparation stage in order to squeeze the right quadra- 
ture of the vacuum. The time delay of the feedback path 
must be much shorter than the time scale at which the 
phase modulation varies. 

If the feedback of 0' and the ^'-quadrature squeezing is 
not feasible in practice, another option is to use an ideal 
phase-squeezed state and squeeze the yg quadrature of 
the vacuum before displacement and phase modulation. 
The squeezed- vacuum Wigner distribution is 



Wo[xo,yQ] cxexp 



-xo-K^^ -xo- -yo -K^^ -yo 

(4.63) 
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FIG. 5: (Color online). A homodyne phase- locked loop that 
realizes the MAP estimation for squeezed states. S{z') de- 
notes squeezing of the z' quadrature. 



The MAP equation becomes 



m = Km H r], 



(4.64) 



where 



rj = q{K-' ■ p) + p{K-' ■{2\a\- q)) , (4.65) 

q = aexp{—i4>) + a* cxp{i(p), (4.66) 

p = — ?'[aexp(— i0) — a* exp(i0)]. (4.67) 

The MAP equation now depends quadratically on two 
quadrature fields q and p. Realizing this MAP esti- 
mation would require some nontrivial second-order field 
measurements of r/ in a Weyl-ordered operator form for 
the statistics to obey the Wigner distribution. 

If we insist on using the homodyne phase-locked loop 
for phase-squeezed states, the noise term of the homo- 
dyne detection z' in the first order of small {{(f) — </>')^) 
is given by 



z' W Xo{(j)-4>') 



Vo- 



(4.68) 



The ijq quadrature is squeezed, but the antisqueezed Xq 
quadrature also enters the phase-locked loop through the 
phase-sensitive term xo{(p— cf) ), the magnitude of which 
varies depending on the difference between the mean 
phase and the local-oscillator phase. If the magnitude 
of Xq{(P — (p ) is much smaller than the magnitude of t/g, 
however, we can assume that z' is approximately j/q, 



Vo, 



(4.69) 



and the homodyne phase-locked loop can still be used, as 
shown in Fig. [SI with the same filters as those designed 
for the G given by Eq. (gUl). 



Because z' passes through the loop filter L before en- 
tering the system, the approximation z' « j/q requires 



((</'-0r)« 






(4.70) 



If the constraint given by Eq. (|4.70p is violated, z' can- 
not be regarded as a phase- insensitive noise term, and the 
phase-locked loop is no longer an accurate realization of 
MAP estimation for phase-squeezed states. It might still 
be possible to heuristically adjust the parameters of the 
phase-locked loop to compensate for the phase-sensitive 
noise and minimize the estimation error, but Eq. (|4.70p 
sets a rough limit on the right amount of squeezing, be- 
yond which stronger squeezing is no longer useful. 

If one is able to perform canonical temporal-phase 
measurements, then the homodyne detection in the 
phase-locked loop shown in Fig. [S] can be replaced with 
measurements of sm{<p — </>), where (p is the canonical 
phase of the incoming field a. Canonical phase mea- 
surements are not sensitive to the antisqueezed photon- 
number noise, so the squeezing constraint given by 
Eq. (I4.70p is not needed, although the constraint given 
by Eq. (|4.4ip is still required for the estimation to remain 
efficient in the linear regime. 

The adaptive temporal-phase measurement scheme 
proposed by Berry and Wiseman [2^ is essentially a first- 
order homodyne phase-locked loop, with an integrator as 
the loop filter and no post-loop filter. While the first- 
order phase-locked loop is indeed the optimal demodula- 
tor when the mean phase is a Wiener random process [4] , 
it is clear from the preceding discussion that, at least for 
stationary messages, one should use a more complicated 
loop filter as well as a post-loop filter according to the a 
priori message statistics to minimize the estimation er- 
ror. In this paper we focus on stationary messages and 
the use of MAP estimation to design angle demodula- 
tors. For nonstationary messages, it is more appropriate 
to use a state- variable approach and Kalman-Bucy filters 
[1, [33I , although it is beyond the scope of this paper to 
study the latter case. 



V. QUANTUM LIMITS TO ANGLE 
DEMODULATION 



To calculate the mean-square error in the estimated 
message by MAP estimation, consider the linearized 
MAP equation 
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FIG. 6: (Color online). A homodyne phase-locked loop for 
phase-squeezed states. 



rh = G 4> = G H m + G —-. (5.1) 

2\a\ 

The mean-square error is 
((m - mf) « Diag \{G-H-I)- K^ {G H if 



4|a| 



^G K2 G' L 



(5.2) 



where Diag(A)j = Ajj is the diagonal vector component 
of A. Equation (|5.2p can be further simphfied in the 
Fourier domain to give 

a^ = {{fhj - rrijY) 

"/ .i„i9o ,M^^^T,M^^■) . o ^jn ' (5-3) 



B 



-B/2 



A\a\^Srr.{f)\H{f)\^+S,{f)' 



where S2{f) is the power spectral density of the z' 
quadrature. This expression is well known in linear es- 
timation theory and called the irreducible error [l|, [3]- 
Assuming that the message has unit variance (m^) = 1 
and a flat spectral density with bandwidth 6, 



Sr 



.if) = { 



B/b, I/I < 6/2, 
0, I/I > 6/2. 



(5.4) 



If S2{f) is also flat over the bandwidth b, the mean-square 
error for PM becomes 



"■pM 



i> /?^A»1, (5.5) 



/?2A+1 /?2A' 



where we have defined the parameter A as 
4|a|25^(0) 



A; 



5*2(0) 



(5.6) 



The SNR is then 



SNRpM = ^^ « (3^A. 






(5.7) 



PM 



For FM, defining the equivalent modulation index as 



P- 



b ' 



(5.8) 



the mean-square error is 
o-pM = 1 - \/^52Atan"i —= 



and the SNR is 



y^ 3/?2A 



SNRfm ~ 3/3^A. 



, I3^K > 1, 



(5.9) 



(5.10) 



The parameter A depends on the ratio of the message 
spectral density to the quadrature-noise spectral density. 



A. Standard and Heisenberg quantum limits 

For a coherent state with average optical power 

V = hhB{a)a) = hfoB\a\^ (5.11) 

and a quadrature-noise spectral density 

S2{f) = 1, (5.12) 



A is 



4P 
A = ^ = 4AA, 



M = 



V 



(5.13) 



hfob ' hfab' 

where M is the average number of photons in a period of 
1/6, and the ct^ ~ 1/N dependence is an analog of the 
usual standard quantum limit (SQL) in single-parameter 
estimation 10]. We can then write the SNR at SQL as 

SNRpM ~ ^SNRfm ~ ^P^M. (SQL) (5.14) 

This linear dependence on A/" may not hold for nonban- 
dlimited message spectral densities. For example, if we 
let the message spectral density be Lorentzian, 

^rr^if) = ^/2 + (5/2)2' (5-15) 

the PM SNR for a coherent state in the limit of _B ^ 6 



SNR 



PM 



8/3W 



1, 



(5.16) 



which scales with VA instead of TV. A similar VA scal- 
ing is also observed by Berry and Wiseman for their adap- 
tive temporal-phase measurement scheme with the non- 
stationary Wiener process as the mean phase [2^. In 
the following we consider only the bandlimited message 
spectral density given by Eq. (j5.4p for simplicity. 
For squeezed states, it is shown in Appendix [Bl that 



52(/)=exp(-2r)for|/|< 



B. 



(5.17) 



where r is the squeeze parameter and Bg is the squeeze 
bandwidth. We should therefore make Bg at least as 
large as the message bandwidth 6. The average power of 
a squeezed state is 



V = hfQB\a\^ + hfoBs sinh^ r. 


(5.18) 


ssuming Bg — b, 




A = 4(7V - sinh^ r) exp(2r). 


(5.19) 


e SNRs become 





SNRpM ~ -SNRfm ~ ^0^(M - sinh^ r) exp(2r). 



(5.20) 



The optimal exp(2r) at which the SNRs in Eq. ((5?20|) are 
maximum is given by 

exp(2r) = 27V+1, (5.21) 

and the maximum SNRs become 

SNRpM « ^ SNRfm « 4/3W(AA +1). (Heisenberg) 

(5.22) 

The N"^ scaling can be regarded as an analog of the 
Heisenberg limit in quantum single-parameter estimation 



m- 
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B. Threshold constraint 

It must be stressed that the preceding derivation of 
quantum hmits is accurate only in the Unear regime, 
when the constraint given by Eq. (|4.4ip is satisfied. Be- 
yond the hnear regime, the phase-locked-loop demodula- 
tor is no longer an efficient estimator, and the SNR is ex- 
pected to drop abruptly when the quantum noise is larger 
than a certain threshold. Physically, the abrupt drop in 
SNR below threshold is due to the periodic nature of p' 
with respect to (f>^ (f> , which can cause to differ from 
(p by multiples of 2tt when the noise is large. To ensure 
that the estimation is efhcient, Eq. (|4.4ip . which we now 
call the threshold constraint in accordance with classi- 
cal theory, should be observed. Since we have assumed 
that the ^'-quadrature spectral density is flat over the 
message bandwidth, an expression for {{(f>j — 01)^) can 
be borrowed from well known results for classical signals 
with additive white Gaussian noise [l|, 0- For the flat 
message spectral density. 



{{^,-^',r) 



1 



6/2 



-b/2 



df\n[l + A\H{f)\' 



The threshold constraint for PM becomes 



-In(l + /32A)<1, 
while the constraint for FM is 



(5.23) 



(5.24) 



ln(H-/32A) + 2/3VAtan- 



PVA 



< 1. (5.25) 



In practice, it has been determined numerically for clas- 
sical signals with additive white Gaussian noise that 
CTn < 0.25 is sufficient to satisfy the threshold constraint 
[l|, [4I , so a similar constraint on (Tq can be used for co- 
herent states or 2;'-quadrature-squeezed states. 



for FM, where A is given by Eq. (|5.19p . The squeez- 
ing constraint is much more stringent than the threshold 
constraint and more so for stronger squeezing, because 
the former constraint requires the enhanced noise in the 
antisqueezed Xq quadrature to be negligible. 

The squeezing constraint also prevents the homodyne 
phase-locked loop with phase-squeezed states from reach- 
ing the Heisenberg limit given by Eq. (|5.22p . Con- 
sider, for example, the left-hand side of Eq. (I5.26P at 
the Heisenberg limit. 



(2AA+1)^ 



ln(l + 4/?2A/-2), 



(5.28) 



which is much smaller than 1 only when AP'^M'^ ^ 1 and 
the exact SNR approaches its lowest possible value 1, the 
a priori SNR. If we let exp(2r) be a small parameter A 
times Af instead. 



exp(2r) = W, 
Eq. (|5.26p becomes 

A< 



A<1, Aw4A7V2, (5.29) 



(5.30) 



ln(l -I- 4A/3W2 



In A/" 



for 4(5/3^2 > 1 and In TV > ln(2^/A/3). The SNR for 
PM is therefore limited by 



SNRpM ~ ^Xp'^Af^ < 



8/3W 
In A/" 



(5.31) 



in the limit of large M . The limit on the FM SNR can 
be derived using a similar argument, and is given by 



SNRfm < 



24/32 AA^ 
In A/" 



(5.32) 



in the limit of large N. The Af^/hiAf scaling is analo- 
gous to the ultimate limit to adaptive single-mode phase 
measurements urn. 



C. Squeezing constraint 



D. SNR enhancement by wideband angle 
modulation 



For phase-squeezed states and the homodyne phase- 
locked loop depicted in Fig. [6l one should also apply the 
squeezing constraint (|4.70p for the SNR to be given by 
Eq. (|5.20p . The order of magnitude of the right-hand side 
of Eq. (|4.70p is calculated in Appendix [B] and is roughly 
exp(— 4r). The constraint becomes 



exp(4r) 
A 



ln(l + P^A) < 1 



for PM and 

exp(4?') 
A 



ln(l + P^A) + 2/3VAtan-i [ ^^ 

V/3VA 



(5.26) 

<1 

(5.27) 



Apart from the quantum enhancement by squeezing, 
the SNRs above threshold can be enhanced quadratically 
for a bandlimited message spectral density simply by in- 
creasing the modulation index /3. This SNR enhancement 
is achieved at the expense of bandwidth resource, as is 
well known in classical communication [l|, l^j ISi- The 
SNR enhancement effect for discrete FM in the quan- 
tum regime has also been suggested by Yuen [3| . A good 
approximation of the optical signal bandwidth B when 
the phase or instantaneous frequency of a monochro- 
matic wave is modulated by a message with unit vari- 
ance (m3) — 1, bandwidth b, and modulation index B is 

1—1 

provided by Carson's rule [3| , 



B«(/3+l)6. 



(5.33) 
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The modulation is described as narrowband when /? <C 1, 
wideband when /3 '^ 1, and ultrawideband when /3 3> 1. 
For squeezed states, because they aheady have an in- 
herent bandwidth Bg before phase modulation, Carson's 
rule should be modified to read 



B^Bs + iP+l)b. 



(5.34) 



For a given A, /3 is also limited by the threshold con- 
straint given by Eq. (jgTM)) . ((05)) . dElSD, or (jOT)) . since 
CTq increases logarithmically with /3. Thus /? cannot be 
increased indefinitely even if enough optical bandwidth 
is available. 



VI. 



OPTICAL SENSING OF POSITION AND 
VELOCITY 




FIG. 7: (Color online). Multipass position and velocity sen- 
sor. 

The preceding analysis can be applied directly to opti- 
cal sensing of position and velocity, if the phase modula- 
tion is due to reflection off a classical moving target. Con- 
sider the setup shown in Fig.[71 A light beam, which can 
be a coherent state or squeezed states, shines on a reflect- 
ing target and bounces back and forth between the target 
and a perfectly reflecting mirror for M times before leav- 
ing the apparatus. The phase modulation due to the 
movement of the target is then estimated using canoni- 
cal temporal-phase measurements or a phase-locked loop. 
The mean phase of the output light beam is given by 



(/){t) = (2Mcos0)^x{t) 



(2M cos ( 



27r/o 



dt'v{t'). 



(6.1) 
(6.2) 



The quantum limits derived in the preceding section are 
also applicable here, if we assume that the target position 
and velocity are stationary random processes and define 



the normalized PM and FM parameters in terms of the 
sensing parameters as 



m{t) 



(3 = 2M cos e 
x{t) 



2Tr^^H^ 



for position sensing, and 



T^{2Mcose) 



f3^ {2McosB) 



c 
be 



l{t) = 



,{t) 



VWi^ 



(6.3) 
(6.4) 

(6.5) 
(6.6) 
(6.7) 



for velocity sensing. The power spectral densities of the 
position and velocity are related to each other by 



S.if) = 



1 



(2^/)^ 



r^.(/), 



(6.8) 



so one should calculate the mean-square errors for posi- 
tion and velocity estimations separately using Eq. (|5.3|1 . 
It is apparent from Eq. (|5.3p that increasing (3 and there- 
fore H{f) can at best reduce the mean-square errors of 
the estimations quadratically Even though the parame- 
ters (a;^(t)), (w^(t)), and h are given for each target, it is 
still possible to enhance the SNR by increasing M, the 
number of times the target is interrogated. In practice 
the maximum achievable M is limited by the threshold 
constraint as well as other experimental constraints. For 
instance, the total interrogation time of the light beam 
with the target should be much shorter than the time 
scale at which the position and the velocity changes, or 
in other words, 



2(Af - l)L 1 

— ^ — < -. 

c cos 6 b 



(6.9) 



The target size, shape, and reflectivity also limit the max- 
imum achievable M in practice. 

The multipass setup may be regarded as a continuous- 
parameter generalization of multipass single-parameter 
quantum estimation schemes [10, [ll|. For single- 
parameter estimation, the multipass SNR enhancement is 
achieved at the expense of time, whereas for continuous- 
parameter estimation the enhancement effect utilizes 
both time and bandwidth resources. Squeezing becomes 
useful when such resources are limited, and the experi- 
ment can be highly controlled to eliminate decoherence, 
such as that caused by the imperfect reflectivity of the 
target or other optical losses in the system. 

To achieve a large effective M, one may also use a 
Fabry-Perot setup, as shown in Fig. [51 similar to the 
interferometric gravitational wave detector with Fabry- 
Perot arms i3!l. cos 6* becomes 1, and the effective M at 
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FIG. 8: (Color online). Fabry-Perot position and velocity 
sensor. 



resonance is 



M = 



R 



1-Vr' 



(6.10) 



Because the resonance condition is sensitive to the phase 
shift by the target, the Fabry-Perot configuration is use- 
ful only in the narrowband modulation regime /3 ^ 1. 
Apart from this restriction, the preceding discussion on 
the multipass configuration applies to the Fabry-Perot 
setup as well. 



VII. FLUID VELOCITY OPERATOR FOR 
NONRELATIVISTIC BOSONS 

We now apply the time-domain optical formalism de- 
scribed in Sec.HIland lllll to nonrelativistic bosons in three 
spatial dimensions. In momentum space, the boson an- 
nihilation and creation operators obey the commutation 
relation 



[a{k),aHk')]^6^k-k'). 



(7.1) 



If we regard the nonrelativistic theory as an effective the- 
ory that ceases to be accurate beyond a certain range of 
momenta, we can explicitly impose a momentum cutoff 
to the Hilbert space, analogous to the bandwidth limita- 
tion on optical fields in Sec.HIl 



l^a;|,|fcyUfcz| <K, 

SO that we can discretize space as 



(7.2) 



Xj=xo+j5a;, i^je^ + key + le^, Sx=—. (7.3) 

K 

We can then define the space-domain operator as 



i(x) = 



1 



K 



dkr 



K 



X a(k) exp(ik • x) 






K 



dk. 



-K 



(7.4) 



and the discrete wave-packet-mode annihilation operator 
as 



aj^i(xj)(fe)3/2, 
with the commutator 



[aj,4] 



-^jj'- 



The Hilbert space can be spanned by Fock states. 



E 



%! 



■|0)j 



(7.5) 
(7.6) 

(7.7) 



Unlike photons, the total number of massive bosons is 
conserved and usually fixed, so it is physically justifiable 
to impose upper limits to the boson numbers, so that the 
Hilbert space consists only of finite-boson-number states. 



Ei")(^ 



(7.8) 



where s = {. . . , Sj, . . . } and for simplicity we let Sj be 
the total number of bosons N . 

With the finite-momentum and finite-boson-number 
Hilbert space, we can now apply the Pegg-Barnett theory 
p^ and define the unitary exponential-phase operator as 



N 



exp(i(^j) = ^ |nj - l)j(nj| + exp[i(A^ + l)0oj]|A^)j 



n-,=0 



(7.9) 



Most importantly, the action of this operator on the vac- 
uum state is uniquely defined as 



exp(z(^j)|0)j=exp[z(A^ + l)(/.oj]|7V)j, 



(7.10) 



which makes the operator unitary and distinguishes it 
from the nonunitary Susskind-Glogower operator. 

The creation and annihilation operators can be rewrit- 
ten as 

flj = exp(i0j)y^, a] = y^exp(-i0j). (7.11) 

Since the Hamiltonian is a function of fij and aj, the 
unphysical jump from the vacuum state to the maximum- 
number state indicated by Eq. (|7.10p due to the action of 
the Pegg-Barnett operator cannot occur directly in the 
dynamics. 

The commutation relation between the Pegg-Barnett 
operator and the number operator fij = " ' 

[exp(z^j),rij] = [l-(N+l)\N);{N\]eM^h)- (7-12) 

To define a fluid velocity operator, consider the semiclas- 
sical definition of a superfluid velocity field [23| , 



a-a\ IS 



V(x) EE -V, 

ra 



(x), 



(7.13) 
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where m is the mass of each boson. Because (/'(x) can 
be a multivalued function and is undefined where the 
fluid density is zero, Eq. (|7.13p is ill-defined even in the 
semiclassical regime. An alternative is to consult the 
definition of the current density [23, [2a] , 



we have 

[v(x), /5(x')] « -i^V^S^i^ ~ x') cos[0(x') - <^(x)] 

«_iAv,j3(x-x'), (7.25) 



J(x) = -- v4*(x)Vxv4(x) - A(x)Vx^*(x) 
Aim L 



(7.14) 



and define the analogous fluid velocity in terms of the 
exponential-phase function, 



v(x) = 2— |exp[-i0(x)]Vxexp[i0(x)] 
- exp[i(/)(x)]Vx exp[-i(/)(x)]| 



-Vx' sinf(/)(x' 



(x)] 



(7.15) 
(7.16) 



which is equivalent to Eq. (|7.13p where (?!)(x) is continu- 
ous, and the sine function is always single-valued where 
(/)(x) is defined. But we have not yet solved the prob- 
lem of defining 0(x) where the density is zero. This can 
be done by regarding 0(x) as a random function, or in 
the quantum regime, an operator. In the discrete space 
domain, we can define the fluid velocity operator as 



j' 



i(0j/ - <^j), 



(7.17) 



the last expression of which agrees with the commutation 
relation proposed by Landau [2J|. Thus, we have shown 
that it is possible to rigorously define Landau's fluid ve- 
locity operator, if we start with the physically reasonable 
assumptions of finite momentum and finite boson number 
and take the continuous limit at the end of a calculation. 
If the continuous limit is taken, one does not have to 
use the Dj_j' defined in Eq. (|7.19p . and any Dj_j' that 
possesses the continuous limit given by Eq. (|7.23p will 
suffice. 

The fluid velocity is a physical quantity that can be 
measured optically. The bosonic fluid can act as a moving 
dielectric, in which a propagating light beam can acquire 
a phase shift proportional to the fluid velocity in the first 
order, due to the Fresnel drag effect 40]. Bosons may also 
reflect or scatter light, and the fluid velocity will shift the 
instantaneous frequency of the reflected or scattered light 
due to the Doppler effect. Thus, the fluid velocity and 
the optical instantaneous frequency are closely related 
physical quantities, and, as we have shown in this paper, 
can be described by the same theoretical formalism in 
the quantum regime. 



where 

sin((^j' 



-^j) 



1 r 

2i 



exp(i(? 



j)-H.c. 



(7.18) 



Dj_j/ is the discrete impulse response that corresponds 
to the gradient operator, 



Di 



1 

bx 



{dj-j'Skk'Sii'ex + dk-k'Sjj'Swey 



(7.19) 



where dj-j' is the differentiator given by Eq. (I3.13p . 
Neglecting the ]A^)j(A^] term in the commutator in 
Eq. (j7.12p , which is important only in the highly unlikely 
event that all bosons arc in the same spatial wave-packet 
mode, we arrive at the following velocity-number com- 
mutation relation, 



n. 



[vj,ny] « -i—D._y cos((/)j' - 0j). 



(7.20) 



In the limit of i5a; ^ 0, 



cos 



Vj -* v(x). 


(7.21) 


nj — > Sx p{x.), 


(7.22) 


Dj_j, ^(5a;3Vx5^(x-x'), 


(7.23) 


-0j)^cos[0(x')-0(x)], 


(7.24) 



VIII. CONCLUSION 

In conclusion, we have proposed a quantum theory of 
optical phase and instantaneous frequency in the time 
domain. In the formalism we have introduced an ex- 
plicit optical bandwidth cutoff, in order to reflect our 
experimental limitations, satisfy the slowly-varying en- 
velope approximation, and, most importantly, discretize 
time domain into discrete modes, the phases of which we 
know how to define and measure in principle. Guided by 
insights from classical estimation theory, we have sug- 
gested the use of homodyne phase-locked loops to per- 
form angle demodulation, and the quantum limits are 
derived. We have also shown how the SNR enhancement 
effect of wideband angle modulation can be applied to 
optical sensing of position and velocity. Given the re- 
cent experimental advances in quantum-enhanced mea- 
surements [lO, [U SH and the maturity of optical phase- 
locked loop technology [42], we expect our theoretical 
predictions to be experimentally realizable with current 
technology and relevant to future communication and 
sensing applications. Finally, we have applied the optical 
formalism to nonrelativistic bosons and shown how Lan- 
dau's fluid velocity operator can be deflned rigorously, 
thus resolving a long-standing issue in quantum hydro- 
dynamics. A more rigorous formulation of the current- 
algebra approach to quantum hydrodynamics based on 
our theory can be envisaged. 
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APPENDIX A: ALGEBRA IN THE DISCRETE 
TIME DOMAIN 

For convenience we adopt the following simplified no- 
tations to describe algebra in the discrete time domain. 
Boldface lowercase letters denote vectors in the discrete 
time domain, 



a = {. . . ,aj,aj+i, . . .}. 



(Al) 



Multiplication and division of components of two vectors 
are written implicitly, 



ab^ {..., ajbj,aj+ibj+i, . . . }, 



a 
b 



' bj' bj+i 



(A2) 
(A3) 



A vector of functions of each component of a is written 
as a function of the vector, 

/(a)^{...,/(a,),/(a,+i),...}, (A4) 

while a function of all components is written with square 
brackets, 

/[a] = /(..., aj,aj+i,...). (A5) 

Dot product of two vectors is defined as 

a-b="^ajbj. (A6) 



The tensor product of two vectors is 
{a<g) b)jk = a-jbk- 



A matrix Ajk is written as boldface and uppercase A. 
The transpose of a matrix is defined as 

{A^),k^Ak,. (A8) 

Dot products of a matrix and a vector are 

{A-b)j=J2^i'^bk, (b.A)j = ^6feAfc,. (A9) 



Multiplication of two matrices is 



Ik- 



(AlO) 



The inverse of a matrix, if it exists, is A , 

A A^ = A^ A = I, (All) 

where Ijk = Sjk is the identity matrix. 



APPENDIX B: BROADBAND SQUEEZED 
VACUUM 



We shall study the statistics of the squeezed vacuum 
|19l |43| in the interaction picture. For simplicity we as- 
sume that the squeezing spectrum is uniform and has a 
bandwidth of Bs- Starting with the vacuum quantum 
state, the output annihilation operator in the frequency 
domain of a bandlimited squeezed vacuum can be written 
as 



i(/) = { 



Mavac(/)+i^atac(-/), \f\<Bs/2. 



avac(/), 



\f\>Bs/2, 



(Bl) 



where fi and v are parametric gain parameters that sat- 
isfy 



IMI- 



(B2) 



and avac(/) is the annihilation operator with vacuum 
state statistics. Performing the inverse Fourier transform 
on Eq. dHI]), 

a, = ^( dfb{f)eM-^27Tft,), (B3) 

VB J-B/2 



1 f^/^ 

flvacj = ^= / 4favac(/)cxp(-l27r/ij), 

VB J-B/2 



(B4) 



the discrete-timc-domain operator for a squeezed vacuum 



(^ - l)avac + i^ata 



where Tjk = Tj-k is a low-pass filter defined as 



Tj-k = 



B 



Bs/2 



Bs/2 



dfexp[~i2TTf{tj ~tk)]. 



(A7) We define the quadrature operators as 



a^o = a + a' , 



Vo 



-i{a — d' 



(B5) 



(B6) 



(B7) 



The vacuum-state operators have the following statistics. 



(ala 



\^vac/ 
a vac) 



= 0, (dv 



0, 



alac) = I, 



which lead to the following covariance matrices, 

Ki = {xo ®Xo)=I-{l-\ti + v*\^ )r, 
K^^ {y^®y,;) =1 - {l-\^i-v*\')T. 



(B8) 
(B9) 

(BIO) 
(Bll) 



To produce squeezing in the ^q quadrature, both /i and 
v should be real and positive. Writing /i and v in terms 
of the squeeze parameter r. 



jjL = cosh r. 



V — sinh r, 



(B12) 



The covariance matrices become 



Ki = / - [1 - exp(2r)] r, (B13) 

K2= I - [l- exp(-2r)] T. (B14) 



The quadrature power spectral densities are 

, , r exp(2r), I/I <Bs/2, 
Si{f)~[ ^^ |/|>S./2, ^^^^' 

, , r exp(-2r), I/I < B,/2, 
S2{f)-[ , \f\>BJ2. ^^^^' 



If Xq and yo P^^ss through a low-pass filter L with band- 
width b < Bs, their variances become 
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and their ratio is 



((i-ioH 



exp(— 4r), 



(B19) 



((L-Ao)'> = -exp(2r), 
((i-yo)2) = Acxp(-2r) 



(B17) 
(B18) 



which provides an estimate of the magnitude of the right- 
hand side of the threshold constraint given by Eq. (|4.70p . 
if we regard i as a bandpass filter with bandwidth h. 
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